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CN ' Abstract 

In a previous paper we treated within the framework of our Projective Unified Field Theory (Schmutzer 
2004, Schmutzer 2005a) the 2-body system (e.g. earth-moon system) with a rotating central body in a 
rather abstract manner. Here a concrete model of the transfer of angular momentum from the rotating 
central body to the orbital motion of the whole 2-body system is presented, where particularly the transfer 
is caused by the inhomogeneous gravitational force of the moon acting on the oceanic waters of the earth, 
being modeled by a spherical shell around the solid earth. The theory is numerically tested. 

Key WORDS: transfer of angular momentum from earth to moon, action of the gravitational force of the 
^ ' moon on the waters of the earth. 
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^ ' 1 Introduction 
in 

■ Nowadays the concept that the tidal braking effect of the earth rotation is caused by the gravitational force of the 
moon acting on the viscous oceanic waters around the earth is accepted in general. With respect to the details 
of this braking mechanism several models were proposed and theoretically investigated in detail, including 
numerical computer calculation. We refrain from reviewing the historically first ideas and the afterwards 
following publications on this subject, but we rather cite some previous papers, where the interested reader can 
find further information: Brosche 1975, 1979, 1989 as well as Brosche and Schuh 1998. 

Our considerations start with a strongly simplified model of the earth consisting of a rigid sphere surrounded 
by a homogeneous viscous fluid shell (oceanic waters) of uniform thickness, being under the influence of the 
rotation and the gravitational field of the rigid earth as well as of the gravitational field of the orbiting moon. 
It may be a good help to the reader to have a look on the figures of forming of the tidal bulge mechanism 
without and with phase lag in both the monographs: Lambeck 1980, Sabadini and Vermeersen 2004. For this 
publication the physical situation of the earth-moon system with phase lag is sketched in Fig. ^ 

For the concrete calculations in the following sections we use Fig. [5] (axis of rotation — > z-axis) and Fig. [3] 
(axis of rotation — > perpendicular to the x-y-plane). 

2 Friction force and friction torque 

Now using Fig. 2, we apply Newton's friction force law to the infinitesimal azimuthal component of the friction 
force dFf ric \$ within a viscous fluid (j) viscosity), referring to an infinitesimal areal element dA^ positioned 
at the point with the cylindrical coordinates (g = ro cos a, <S>^ z — ro sin a) at the bottom of the ocean: 

. , _ dv , 2 dv 

a) ah r ric \<f, = n— aA<p = r n — cos a a<P aa with 
1 dr dr (1) 

b) dA$ — g d<P ds , c) ds — ro da , d) v = Log = tquj cos a 

(ro radius of the earth, v radial-dependent azimuthal velocity, to rotational angular velocity). 
In vectorial writing the corresponding infinitesimal friction force reads: 

dFf ric \& — e$ dFf ric \$ . (2) 
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Figure 1: Deformed fluid shell (deformation 
caused by the gravitational field of the moon) in- 
cluding the rotation of the earth, which causes the 
phase lag by the dragging effect within the viscous 
fluid. 




Figure 2: Rigid earth, ocean as viscous fluid shell 
and line element ds at the interface between both 
parts of the earth (tq radius of the rigid earth) . 
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Figure 3: Equilibrium of the tidal bulge (quasi-mechanical body with mass mt in the point Pi) as a result of the 
acting gravitational force of the moon and the dragging force (negative friction force) within the fluid, caused 
by rotation of the earth. 
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Because of the rotation-symmetry assumed we find by integration over the azimuthal angle <P from to 27r the 
intermediate results 

&A$ = 2irr?, cos a da , b) dFr ric = r\— dA . (3) 

dr 

o 

For simplification of the problem we suppose a linear radial velocity profile (d depth of the ocean): 

a) v(r, a) = v{r 0l a) 1 - 1 J oc (r - r ) , i.e. b) = - 1 J° C v(r ,a) , (4) 

where 

v d (a) = s oc v(r ,a) (5) 

is the velocity at the surface of the ocean. We name the constant s oc < 1 "oceanic flow slope parameter" . 
Eliminating the velocity gradient in (PJd), we find 

dF frlc = _ q(l-'~Mro,a) dA (6) 

and further with the help of and (JTJl) 

27rr|o^7(l - s oc ) cos 2 a 

dF fric = da . (7) 

7T 7T . 

Integration over a from — — to — yields 
7T 2 r^??(l - s oe ) 

^/ric -3 ■ (8) 

Adapting this uniform oceanic shell modeling to the earth, where about 70.8 % of the surface are oceanic waters, 
it is convenient to introduce a rough "surface correction factor" f oc . Then formula (JHJ) reads 

_ 7r 2 rg^?7(l - s oc ) f oc ,„s 

" fric i ■ 

As an example one is tempted to choose f oc ~ 0.6. 

Now we define in the usual way the infinitesimal vectorial friction torque 

dM fric \$ = e r r x dF fHc \0 = -e g r dFf ric \$ (10) 

corresponding to the infinitesimal friction force 

In analogy to the above calculations we receive by integration 

f i „ „- , ^ , 2Trriu!n(l — s oc ) cos 3 a da ,„„, 
dM /nc = / dM fncl <p d<P = -k (11) 

and by further integration 

2 

2TTr*UJ7l(l - S oc )foc f 3 , , 1r) s 

Mf r ic = —k / cos a da (12) 

<*=-$ 

or explicitly 

^ A/f un/r m aj 8nr%ujri(l - s oc )/ oc 

a) Mf ri c = KMfric With b) Mfric = ^ ■ (13) 
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Let us now remember some results from our previous paper (Schmutzer 2005a): We used the following definition 
of the rotational angular momentum of the earth rotating about the z-axis: 



L = kL rot = k9 E uj E (14) 

(6 E moment of inertia, u> E rotational angular velocity, k unit vector in z-direction). Hence by temporal differ- 
entiation follows (dot means total time derivative) 

a) = k(6 E uj E + 6 E Lu E ) with b) e E =°Q E + Q E E. (15) 

From this last equation we remember that the total time derivative of the moment of inertia consists of two 
parts: the partial derivative (denoted by circle) and the temporal cosmological influence (determined by the 
"logarithmic scalaric world function" U). 

According to the balance equation of the angular momentum 

^ = M fric (16) 

the temporal change of L is caused by the friction torque (flrtb ,). Inserting the expression 1)15^ and $13h) into 
(f 1 01) yields the following temporal change of the angular velocity (for u> — > lo e ): 

o 

8nr%oj E ri(l - s oc )f oc (9 e + 9 e £)lu e 

" E = seTd T E or (17a) 

- s oc )f oc _ 3[e E dj E + (e E + 9 e S)lo e } 



. i (17b) 

d <AirrQUj E 

3 Deformation of the earth by the inhomogeneous gravitational field 
of the moon 

Basis of our following considerations is Fig. ??. For treating the above physical problem further we first calculate 
the rotation-symmetric gravitational field of the moon (taken as point-like) about the £-axis (R = \/ x 2 + y 2 ) 
under the approximate suppositions: 

a) — -C 1 , b) —2^0 < 1 , c) — ?L° < 1 (18) 
tm msrlj m E r M 

(m E mass of the earth, tum mass of the moon, Tm distance between the centers of earth and moon). 

First we remember that the gravitational force of the moon acting on the center of the earth is compensated 
by the centrifugal force in this same point, caused by the orbital motion of the earth about the common center 
of mass of the 2-body problem. 

Now we begin our considerations on the gravitational potential <f>M of the moon at an arbitrary point P 
inside the earth (r ^ r*o). Using the notation of Fig. ??, then the well-known Newtonian potential reads 

4>m{q) = • 19 

Q 

Hence for the center of the earth results 

4>M(rM) = ■ (20) 

The difference of both the potentials is equal to 

<Pm{q) = <I>m(q) ~ 0M(rjw) = -7n?"m(- —) ■ (21) 



According to the cosine law we learn from Fig. ??: 

1/2 

■ cos<£ 



9 = r M 



{ R\ 2 2R 

1+ cos^ . 22 
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With the help of (|18b ) we approximately find 



Q = r M 
and further 



R _ 1 

cos <P + - 

rn 2 \r M 



7NWM 
TM 



R 



cos<P - -(— Y(l - 3cos 2 <2>) 



Rearranging this equation gives 



a) <Pm(x,v) 



z ' M 



(y 2 - 2x 2 ~ 2xr M ) with b) R 2 = 



y 



(23) 



(24) 



(25) 



With respect to the earth, for simplicity we refer to a homogeneous mass sphere. From usual textbooks we 
take the expression for the gravitational potential in the interior of such a homogeneous spherical body (in this 
notation e.g. Schmutzer 2005b): 



<t>m{r) - 
<t>Ee{r) 



1NTUE f 2 

2rg 
7 N m B 



(r 2 - 3r 3 ) 



(interior potential), 
(exterior potential). 



(26a) 
(26b) 



From $2bh ) and (|26fa ) we receive the superposition potential of the earth and the moon in the interior of the 
earth: 



<t>i {x,y) = &Ei (x,y) + tpM {%, y) 
where the abbreviations 
a) <j) 



lNm 2 M rl 



2m E r% I ( 



1 



m E r M 



2r 
2r 



1 



2m M ro \ (x - x ) 2 



TTl E V' 



M 



A + rn M rl \ y 2 
\ m E r\,l r 2 



b) x 



m M rl 



m E r 2 M (l 



2m M TQ \ 



(27) 



(28) 



m,Er M 



) 



were used. 

Let us here mention that the treatment of this subject in the corresponding literature usually starts from 
the integral representation of the potential. Our superposition method applied above lead us to the resulting 
ellipsoid formula by setting <pi(x, y) = const in (|27[) . This approach shows that under the influence of the moon, 
beyond the deformation of the sphere to an ellipsoid, our calculation exhibits a translation of the resulting 
ellipsoid by the amount xq in the x-direction. 

Using the semi-axes a and b, taken from the above abbreviations, 



a) a 2 



2mur% 



b) b 2 



mMrl 



(29) 



m E r M 
the equation (|2"7Jl reads 

■{x - Xq) 2 



m E r 3 M 



<t>i{x,y) 



7N™£ 



2r 



t-i 

b 2 



2r 



(30) 



Fitting this potential to the surface equation of the physical ellipsoid 



(x - Xq) 2 



b 2 



(31) 



of the slightly deformed earth by the inhomogencous gravitational field of the moon leads to the expression 
7nto e J-Nm E 



tpix^i V) ellipsoid 



2r 



[>0 = 



ro 



2m E r 



M 



1 



2m M rl 
m E r 3 M 



(32) 
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Further simplification of the above formula (|31|) is reached by applying the approximate suppositions derived 
from (JTSJl: 

x m M rl . 37NTOB / m A / rjj \ ~ / m M rj] \ 

a) x = b) <j) = — , c) a = r (H g- , d) 6 = r 1-- §- . (33) 

kieIm 2r V m E r a M ) V 2mMr° M / 

Conservation of the mass of the earth 

we = y Mo^o (34) 
(/io constant mass density) during the deformation leads to the relation 

rl = ab 2 . (35) 

One should remember that for simplicity our above calculations were based on a homogeneous mass model of 
the earth. For our further calculations we are only interested in the physical situation at the surface of the 
sphere, not explicitly referring to its interior. Therefore our further treatment is consistent. 

4 Rotation of the earth, dragging force on the tidal bulges and cal- 
culation of the tidal phase lag angle (tidal dragging angle) 

4.1 Forces on the moon-nearest bulge 

Nowadays there is no doubt among the experts that the phenomenon of the oceanic tides is primarily caused 
by the moon. The observed tidal water bulges in the zenith (direction to the moon) and symmetric to it in the 
nadir (opposite direction) are modeled by the two water bulges resulting from the symmetric deformation of 
the sphere to the ellipsoidal form calculated above. Since the rotational angular velocity of the earth is much 
bigger than the revolution angular velocity of the moon, the earth performs a relatively quick rotation with 
respect to the tidal bulges, i.e. with respect to the moon (motion under the quasi-fixed bulges). Hence follows 
that an observer in the frame of reference of the surface of the earth realizes the phenomenon of tides. 

Here we investigate the dragging effect of the viscous waters, caused in the waters by the rotation of the 
rigid earth. This effect leads to a shift of the bulges by a certain phase lag angle \- The equilibrium position of 
the bulges is the result of two on the bulges acting torques, caused by the gravitational force of the moon and 
the dragging force in the viscous oceanic waters, as above already mentioned. 

Our task is now to determine this equilibrium position of the moon-nearest bulge which abstractively will be 
considered as a mechanical quasi-body, where the forces act on the center of mass of this quasi-body (mt tidally 
caused mass of the quasi-body). The following investigations are based on Fig. ??, where for simplicity the 
ellipsoid is approximated by a sphere and the point of consideration is shifted from the interior to the surface 
of the earth (P -> Pi). 

Let us further mention that the gravitational force of the sun, which is distinctly smaller than the gravita- 
tional force of the moon, will be neglected for simplification of the physical problem, i.e. we concentrate our 
investigation on the primarily essential points of the equilibrium mechanism. 

In order to calculate the resulting torque acting on the center of mass of the moon-nearest bulge, we first 
list the various forces being present: 

F E = m t G E = m t e R G E \ R (36) 
(radial gravitational force of the earth), 

F M = m t G M (37) 
(radial gravitational force of the moon) , 

F c = m t G c = m t eRr Q uj E (38) 
(radial centrifugal force by rotation of the earth), 

Ffric = e<pFf r i c (39) 



G 



(azimuthal friction force determined by @), 

F P = e R F m (40) 

(radial pressure force on the bulge as back-reaction of the earth) . 

One should realize that the quantities Ge, Gm and G c are acceleration quantities. 
Using the radial and azimuthal vectorial decomposition of (|37|l . we obtain 

a) F M = e R F M \ R + e^F M \ $ with b) F M \ R = m t G M \ R , c) F M \$ = m t G M \$ . (41) 

Hence follows the radial component of the total force as the sum of all radial components: 

F R = m t (G E \ R + G m \r + row|) + F P \ R . (42) 

4.2 Stationary equilibrium 

Demanding stationarity we arrive at the equation 

F R = . (43) 

With respect to the azimuthal component of the force we meet the following physical situation: Due to the 
viscosity of the fluid the quasi-body (modeled bulge) is dragged into the direction of the rotational orbital 
velocity by the dragging force 

Fdrag = ~F }nc (F drag > 0, F fnc < 0). (44) 

In the case of stationarity the dragging force is compensated by the (oppositely directed) gravitational force 
(pulling force) of the moon. 

These previous considerations lead us to the equilibrium equation 

F$ = F M \$ + Fdrag = m t G M \$ + Fdrag = (45) 

being by means of (|44|l equivalent to the condition 

Ffric = m t G M \& ■ (46) 
According to the definition of the torque by the azimuthal forces 

M<p = F r , (47) 
for the equilibrium results 

A/<s = , (48) 
being equivalent to (@2J). 

4.3 Calculation of the azimuthal gravitational force component and of the gravi- 
tational torque, both caused by the moon 

After these decomposition procedures we are now able to calculate on the basis of Fig. ?? the azimuthal 
component of the gravitational force of the moon acting on the tidal bulge being at the oceanic surface of the 
earth (R — > Tq, <& — > x)- I n this specialization formula l|23|l reads 



r 1 / rn \ 2 

1 cos ; 



r M i vtm- 

By differentiation follows 

d £ ■ ft , r o 



X+x(— ) (l-cos 2 X )l. (49) 
2 \r M J J 



r sinx(l + — cosx) ■ (50) 
OX V r M i 
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We further receive with the help of l|49|l from l|19|) 

(0 = 1+ cos X 77 I (l-3cos^x) (51) 



'ML. 

r M l r M - v ''.i/ 

and hence by differentiation 



H cosx- I (1-4 cos x) ■ (52) 



dt; r 



By means of (|52f) and (|50|1 follows 

d(f> M d<fr M (£) d£ 7 N m M r sinx 



r 



1 + — cosx- f— ) 2 (1 - 6cos 2 x)j • (53) 



This result leads to the azimuthal gravitational force of the moon, acting on the moon-nearest tidal bulge at 
the azimuthal angle \ (first order approximation in r*o): 

„ r , 9<t>M 7NW M m t sinx/ 3r \ 

F Mx = m t G Mx = - m t tt, r = 2 I H cosx • (54) 

d(r x) r M V r M ) 

Hence by the angle transformation x X + f we arrive at the corresponding force on the diagonally situated 
second tidal bulge: 

„ 7 N m M m t sinx/ 3r \ 

^m( x +tt) = 3 I 1 "— cosx). (55) 

From (|54|) and l|55|) we find for the azimuthal gravitational force of the moon, acting on both tidal bulges: 

^M|az = -^m x + ^m(x+tt) = 3 sm(2x) • (56) 

The corresponding azimuthal gravitational torque, caused by the moon, reads: 

M M \az = F M \ az r . (57) 

4.4 Tidal phase lag angle (dragging angle) 

This quantity follows immediately from the stationary equilibrium condition (|45H . Inserting l|56|) and J3J) into 
this condition gives the intended result (u> — > loe angular velocity of the earth): 

A ■ fo \ 7r2r r M ti W 1 - Soc)foc , s 1 . 7T 2 r^r^w £; ?7(l - s oc )f oc 

a) sm(2x) = -j-^ or b) X = 7; arcsm ^Hr 1 58 

(x tidal phase lag angle or dragging angle). 

For practical reasons we approximate l|58h) for the case \\\ -C — : 

x = Tr 2 rlr 3 M ujEv(l - s oc )f oc 
67 N m Mm t d 

Further it is useful to introduce the meaningful quantity 

a '?(! — s oc)foc /„„ N 
i>Udal = j (60) 

which we name "tidal coefficient" . 

By means of (|59H we receive the following different form of it: 

_ 67 N m M 

& tidal — o 9 -i j (.Drj 

* r r M^E 
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In this context we remember the relation (|17b ). Neglecting the rebound effect, cosmological effect, etc., from 
this relation mentioned we find 

77(1 - S oc )f oc _ 30EUE\tidal ,g 2 - 

Eliminating in l|60|) with the help of this formula, we obtain the following alternative relation for the tidal 
coefficient: 



W E u> E \ t 



& tidal — —7; 4 • ID'jJ 

5 Numerical evaluation 
5.1 Empirical values 

First from the corresponding literature we list some empirical data on the earth-moon system etc., being useful 
for the following numerical calculations (Gauss system of units): 

7 N = 6.673 • 1CT 8 g^cmV 2 (64) 
(Newtonian gravitational constant), 

r = 6.378- 10 8 cm (65) 
(radius of the earth), 

r M = 3.844- 10 10 cm (66) 
(distance earth/moon), 

m E = 5.976- 10 27 g (67) 
(mass of the earth), 

tom = 7.348 ■ 10 25 g (68) 
(mass of the moon). 

lu e = 7.292 • lO^s- 1 (69) 
(angular velocity of the earth), 

^Ma/ = -6.15-l(r 22 s- 2 (70) 
(tidal braking angular acceleration of the earth), 

Hoc ~ 1-025 gem -3 [IERS Conventions 2003] (71) 
(average mass density of the oceanic waters), 

E = C E = 8.039 • 10 44 g cm 2 (72) 

(polar moment of inertia of the earth) . 
Let us here add two remarks: 

1. The values listed above partly serve as data for approximate information. 

2. The following numerical calculations are performed for rough testing our theory with respect to the order 
of magnitude. 
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5.2 Detailed numerical results by using values of the above list 
From (|61|) we determine the tidal coefficient: 

Stwai= 1.77-10- 27 cm- 2 s- 1 , (73) 

whereas from l|63|) we find 

mtxSudai = 4.91 • HTVm-V 1 . (74) 
Using the value J73J), we arrive at 

m 4 X = 2.77-10 18 g. (75) 

As empirically observed, the tide is coming in approximately 25 minutes after the culmination of the moon. 
This fact corresponds to a phase lag angle of 

a) x = 6-25° = 0.11 rad, i.e. b) Bin* = 0.112. (76) 

Inserting the first value into (|75|) , we find for the mass of the tidal bulge the numerical result 

m t = 2.54- 10 19 g. (77) 

Since we here are mainly interested in orders of magnitudes of some physical quantities, we refrain from refine- 
ments of our rough oceanic shell model. For simplicity we therefore choose for the thickness of the shell the 
estimated average value 

d = 3.8-10 5 cm = 3.8km (78) 

to be found in literature. Further for our refinement parameters we take the simple values s oc — 0, f oc = 1. 
Then from (|60[) results 

a) S U dai = 3 or b) r] = S tlda idm t x ■ (79) 

Using the above numerical values (|73|l . (|76fa ). I|77[) and (J75J), we arrive at the following approximate value of the 
viscosity of the oceanic waters 

r) = 1.87- lQ -3 gcm _1 s _1 = 1.87- 10~ 3 Poise (1 Poise = 1 gcurV 1 ) (80) 

which, considering our rough shell model, seems not to be too far away from empirical viscosity values of water, 
e.g. rj water w 10 -2 Poise at 20° C. 

Finally we estimate the seize of the tidal bulge, up to now roughly treated as a quasi-mechanical body. We 
approximate such a body by a spherical cap whose volume is given by the formula (R radius of the cap, h height 
of the cap) 

V cap = ^ (3i? 2 + h 2 ) . (81) 
Rearranging this equation leads us to the following formula for the radius of the cap: 



R=J**--\». (82) 

V niJ oc h 3 

Taking from literature the approximate empirical value of the height of the tide cap h = 53 cm, by means of 
B71JI and 1)77(1 we receive for the radius of the cap the numerical value 

R = 5.46 • 10 8 cm = 5460 km . (83) 

Thankfully I appreciate interesting discussions with Prof. P. Brosche (Observatory Hoher List) and the 
technical help of Prof. A. Gorbatsievich (University of Minsk). 
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